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Abstract. We construct a linear filter optimised for detecting dark-matter halos in weak-lensing data. The filter assumes a mean 
radial profile of the halo shear pattern and modifies that shape by the noise power spectrum. Aiming at separating dark-matter 
halos from spurious peaks caused by large-scale structure lensing, we model the noise as being composed of weak lensing 
by large-scale structures and Poisson noise from random galaxy positions and intrinsic ellipticities. Optimal filtering against 
the noise requires the optimal filter scale to be smaller than typical halo sizes. Although a perfect separation of halos from 
spurious large-scale structure peaks is strictly impossible, we use numerical simulations to demonstrate that our filter produces 
substantially more sensitive, reliable and stable results than the conventionally used aperture-mass statistic. 



1. Introduction 

Are there dark halos in the universe, i.e. dark-matter halos 
of substantial mass which for some reason fail to produce 
light? This question has recently attracted attention by the de- 
tection of significant and apparently massive peaks in weak- 
lensing surveys which are not seen to coiTespond to galaxy 
groups or cluster s in the optical or X-ray bands (see e.g . 
Erben et all I2OO0I: lUmetsu & FutamasS I2OO0I: iMiralles et all 

m 



2002HErben et alil2003l) . 

Dark-matter halos are embedded into large-scale structure. 
Measuring the inhomogeneities of projected mass distributions, 
lensing observations of halos are overlaid by the lensing sig- 
nal of the large-scale structure in front of and behind the ha- 
los. Being approximately a Gaussian random field, lensing by 
large-scale structure adds peaks and troughs to the signal which 
can be mistaken for halos. It may thus be that part of the 
claimed dark-halo detections are in fact peaks in the random 
weak gravitational lensing signal of the large-scale structure. 

Is there a way to separate the two types of lensing signal, 
and thus to improve the reliability of weak-lensing halo detec- 
tions? The principal difficulty in doing so is the unsharp bound- 
ary between halos and the large-scale structure. Nonetheless, 
we shall describe below the construction of a linear matched 
filter which attempts to suppress the weak-lensing signal from 
large-scale structures and yet leaves the signal from sufficiently 
massive halos for detection. 

Filtering in Fourier space is a flexible and powerful tool 
which has many applications for galaxy cluster lensing, like 
the lensing on the Cosmic Microwave Background radia- 



tion, where there is a stro ng superimposition of diff^ erent 
noise contributions (see e.g. I Seliak & Zaldarriagal ( EoOol) and 
^aturi et al. (2004)), and lensing on background sources. It 
could also be used for lensing on background galaxies. 

Our strategy is as follows. The large-scale structure can 
be considered as composed of dark-matter halos of a broad 
and continuous mass range. Yet, there is a characteristic length 
scale at each cosmological epoch, the nonUnear scale, at which 
the variance of the dark-matter density contrast becomes unity. 
This length scale se parates the small-scale regime where the 
dark-matter power spectrum is dominated by the contributions 
of presumably virialised halos from the large-scale regime on 
which the dark-matter density can be considered as a linear su- 
perposition of linearly evolved perturbation modes. In view of 
weak lensing, this suggests the operational definition of large- 
scale structure lensing which is contributed by the linearly 
evolved matter distribution beyond the non-linear scale, and 
halo lensing, which is contributed by non-linear, gravitation- 
ally bound, virialised structures. 

We are aware that this operational definition may be ques- 
tionable for many applications of weak gravitational lensing. 
For instance, it is well known that it is mainly the non-linear 
part of the dark-matter power spectrum which has now been 
measured with spectacular success in the weak-lensing surveys 
of many groups. Yet, we aim at a substantial more modest goal 
here. Being confronted with peak detections in weak-lensing 
data, we wish to know which of them can be attributed to halos, 
and which of them are likely caused by the large-scale struc- 
ture. 
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In this attempt, we adopt the operational definition of 
the distinction between large-scale structure and halo lensing, 
construct a linear matched filter for halo lensing and model 
the lensing signal by intervening large-scale structures as a 
noise component influencing the shape of the matched filter. 
Matched filters for extracting sources of interest embedded 
into substantial background noise were previously proposed 
for many kind of observations, ranging f rom the analysis of 
Cosmic Microwave Background ma ps (Haehnelt & Tegmark, 
| l996t iTegmark & de qiiveira-Cos'tal I1998l: .Sanz et aU .2001 : 
iHerranz et ali l200 2albt ISchaefer et all 20041) to weak lensing 
tomography faennawi & Spergelll2004l) . 

In Sect. 2, we briefly review those aspects of gravitational 
lensing which are needed for the rest of the paper. We describe 
the filter construction in Sect. 3 and contrast it with the more 
conventional aperture mass in Sect. 4. Simulations described 
in Sect. 5 are used in Sect. 6 to show that our suggested fil- 
ter yields promising results. We summarise and conclude in 
Sect. 7. 

2. Weak lensing quantities 

We briefly review in this section those gravitational-lensing 
concepts which we shall use later. For a more de- 
tailed discussion, we refer th e reader to the review by 
[Bartelmann & Sc hneide3 ( l200ll) . We start with an isolated lens 
whose surface-mass density is 1,(6) at the angular position on 
the sky. Its lensing potential is 



4G D\D, 



Dis j 



dV2(6l') \n\0-0'\ 



(1) 



where Di are the usual angular-diameter distances between 
the observer and the lens, the observer and the source, and the 
lens and the source, respectively. The reduced deflection angle 
experienced by a light ray crossing the lens plane at is the 
gradient of the potential. 



The distortion is described by the two components of the shear, 
1 



(7) 



which are combined into the complex shear y = yi + iya- We 
further use the common abbreviation 



(8) 



Outside critical curves, image ellipticities are determined 
by the complex reduced shear 



7(0) 

l-K(0) 



(9) 



In the weak- lensing limit, /<• «: 1, and the reduced shear ap- 
proximates the shear, ^ »; y, to first order. 

Source and image shapes are quantified by the complex el- 
lipticity 



e = 



a — b 
a + b 



,2iO 



(10) 



where a and b are the semi-major and semi-minor axes of an 
ellipse fitting the object's surface-brightness distribution. The 
position angle of the major axis of the ellipse is The expec- 
tation value of the intrinsic source elUpticity is assumed to 
vanish. 

A sufficiently small source with ellipticity is imaged to 
have an ellipticity 



1 + 8*£. ' 



(11) 



where the asterisk denotes complex conjugation. This equation 
illustrates that the lensing distortion is determined by the re- 
duced shear, which is the only lensing quantity directly acces- 
sible through measurements of galaxy ellipticities alone. 



a(0) = V(A(6l) 



(2) 



The image positions for a source located at are given by the 
lens equation 



0^fi- a(0) 



(3) 



For sources much smaller than typical scales on which the 
lens properties vary, the lens mapping can be linearised. The 
deformation of images with respect to the source is then given 
by the Jacobian matrix 



-72 

72 I -K + yi 



Here, k is the convergence 

k(0)=^ = ^((All +'A22) , 



i.e. the surface-mass density scaled by its critical value 



47tG DiD, 



(4) 



(5) 



(6) 



3. Filter construction 

3.1. General description 

We wish to construct a filter allowing the weak-lensing sig- 
nal of halos to be extracted from weak-lensing observations of 
wide fields. We want the filter to suppress, as well as possi- 
ble, the weak-lensing signal of the large-scale structure which 
can otherwise be mistaken as lensing by halos. As described in 
the introduction, we are well aware that this separation cannot 
fully succeed because of the unsharp boundary between lensing 
by large-scale structure and dark-matter halos. According to 
the operational definition of large-scale structure lensing given 
there, we shall model the weak lensing by linearly evolved 
structures as a noise component to be filtered out. 

Consider the weak gravitational lensing signal of a dark- 
matter halo with amplitude A and angular shape t(0). The mea- 
sured signal S (0) is contaminated by some noise N(0), 



S(0)^At(0)+N(0) , 

whose components are to be specified below. 



(12) 
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In our application, the observed signal is provided by the 
measured ellipticities of the background galaxies, which are 
weakly distorted by the foreground halo and the intervening 
large-scale structure. The expectation value of the image ellip- 
ticities at a given position on the sky is the reduced shear. 



(13) 



There are several noise contributions. First, the gravita- 
tional shear is being measured at the random positions of the 
background galaxies. This adds shot noise proportional to the 
number density of the galaxies. Second, these galaxies are in- 
trinsically elliptical, thus the determination of a single galaxy 
ellipticity is a very noisy measurement of the shear only. This 
adds a noise component proportional to the variance of the in- 
trinsic galaxy ellipticity. Finally, there is the noise component 
which we are primarily aiming at, which is caused by weak 
lensing of intervening large-scale structures. This third source 
of noise will be modelled by the linear dark-matter power spec- 
trum, according to our operational separation between halo and 
large-scale structure lensing. 

These noise components are assumed to be random with 
zero mean and isotropic such that their statistical properties 
are independent of the position on the sky. These assump- 
tions are well justified, since in first approximation the back- 
ground galaxies are randomly positioned and oriented, and 
weak lensing by the large-scale structures is well described by 
an isotropic Gaussian random field. The variances of the noise 
components are conveniently described in the Fourier domain, 
where their correlation functions are 



{Nik'TNik)) = (27T)^S(k' - k)PN(k) ; 



(14) 



the hats above symbols denote the Fourier transforms. 

We wish to construct a linear filter ^{0) which, when con- 
volved with the signal, will yield an estimate Aest for the am- 
plitude A at the position 0, 



)T(6»-6/')dV . 



(15) 



We further want the filter to satisfy two constraints. Fkst, 
we require it to be unbiased, i.e. the average error. 



b = <Aest - A) = A 



(16) 



has to vanish. Second, the measurement noise cr, determined by 
the mean-squared deviation of the estimate from its true value, 

a-'' = ((Aest - Af) ^b^- ^ J \nk)f PN(k)cfk , (17) 

has to be minimal. 

For finding a filter satisfying these two conditions, we 
combine them by means of a Lagrangian multiplier A, carry 
out the variation of L = (T^ + Ab with respect to ^' and thus find 
the function ^' minimising L. The solution of this variational 
minimisation is given by 



4'(^) 



1 



(271)2 



PNik) ■ 



d^k 



m 

PN(k) 



This is the usual result showing that the optimal linear filter 
shape is determined by the shape r of the signal, divided by the 
noise power spectrum. The filter is thus made most sensitive 
for those spatial frequencies where the signal f is large and the 
noise P^ik) is small. 

3.2. Modelling the signal 

Before we can apply this formalism to the observed data, 
we need to assume a model for the spatial distribution of 
the signal and a power spectrum describing the noise proper- 
ties of the data. Concerning the signal, we assume that clus- 
ters are on average axially symmetric. Thus, t(0) = t(\0\). 
Specifically, we use the density pr ofile found in numerical sim- 
ulations by iNavarro et alJ ( 1 19971) for modelling the average 
three-dimensional density profile of dark-matter halos. 



p(r) = 



(r/rs)(l + r/rs)2 



(19) 



where ps and rs are characteristic density and distance scales, 
respectively. These two parameters are not independent, but 
stochastically related to a single parameter, which can be cho- 
sen as the cluster mass. 

The gravitational lensing properties of such a mass distri- 
bution have been widely explored (see e . g . ' B artelmannl [l 9961 
Meneghetti et aL .20031) . Its lensing potential is given by 



tf/(x) — 4Ksh(x) , 
with 

h(x) =-\n^- + 
'2 2 



(20) 




(x>l) 

{x<\) 
(x= 1) 



(21) 



where x is the projected distance from the lens centre in units 
of the distance scale, x - rjr^. The reduced shear g is obtained 
from Eq. (I20> using Eqs. (|5} and 0. Its Fourier transform is 
defined as 



f(^) = ^{k) - I d^x g{x) exp(ijr ■ k) 



(22) 



3.3. A/o/se power-spectrum 

As discussed earlier, we take three sources of noise into ac- 
count. These are the weak-lensing signal due to the large-scale 
structure, the shot noise from the finite number of randomly 
placed galaxies, and the noise due to the random intrinsic shape 
and orientation of the sources. 

The weak gravitational lensing effect of the large-scale 
structures between the observer and the sources at a comoving 
distance w, observed in the directio n on the sky, gives rise 
to an effective convergence (see e.g. lBartelmann & Schneiden 
12001 1) 



(18) K,g(0,w)^ 



2c2 



fd.'^ 

Jo 



(w-w')Mw') 6[Mw')0,w'\ 



a{w') 



(23) 
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The factor out front contains the present-day matter-density 
parameter and the Hubble constant Hq. The comoving 
angular-diameter distance 

C R-^l'^ s\n(K^l^w) K>Q 
Mw) = I w K^O (24) 

is parametrised by the spatial curvature K. The density contrast 
6 = (p - p)/p is taken at perpendicular distance fK(w')0 and 
parallel distance w', and a{w') is the scale factor of the universe 
at distance w'. 

When the sources are distributed in comoving distance, 
Kgff(0, w) needs to be averaged over the source-distance distri- 
bution G{w), 

K,ffi0)^ dwG(w)K,ffie,w) , (25) 
Jo 

where the upper integration Umit wh is the comoving distance 
to the horizon. 

The power spectrum of the effective convergence P^ik) is 
related to the power spectrum of the three-dimensional density 
fluctuations Psik) by Limber's equation, 

4c^ Jo a (vv) XMw) I 

in which the weight function W(w) is given by a line-of-sight 
integral over distance ratios, 

r"" Mw' - w) 

Jw MW) 

The cosmic-shear and effective-convergence power spectra are 
identical, 

Py(k) = PJk) . (28) 

The noise contributions from random galaxy positions and 
eUipticities are both modelled as Poisson noise, which has a 
flat power spectrum. The shot-noise power spectrum from the 
galaxy positions is given by the number density of suitable 
galaxies, 

Pp(k) cc - . (29) 
"g 

The noise component caused by the random intrinsic shape and 
orientation of the background galaxies has a power spectrum 
whose amplitude depends on the variance cr- of the intrinsic 
eUipticity distribution, 

Pe, <x crl . (30) 

The two noise components are combined into one power spec- 
trum which quantifies the noise of an ellipticity measurement, 

Peik) = ^ . (31) 

«g 

With the remaining noise components being independent, the 
total noise power spectrum is the sum of the contributions, 

PN{k) = Py{k) + P,(k) . (32) 




Fig. 1. Power spectra of the various components contributing 
to the definition of the optimal weak-lensing halo filter: effec- 
tive convergence (dashed-dotted line), white noise composed of 
shot noise from the random galaxy positions and their intrinsic 
eUipticities (dotted line), and the model signal (solid line). The 
resulting filter is shown by the heavy solid line. 

The contributions affecting the construction of the opti- 
mal filter are shown in Fig.^ The solid line shows the power 
spectra of the signal, assuming an NFW lens model of M = 
5x10'"^ Mq at redshift z\ - 0.3 and a source redshift Zs - 2. 
The effective shear power spectrum and the Poisson noise are 
given by the dash-dotted and the dotted lines, respectively, 
where the latter contains the contributions from the random 
galaxy positions and their random eUipticities. The bold solid 
line shows the filter resulting from Eq. ( I18> . Unfortunately, 
but expectedly, the large-scale structure power is significant on 
those spatial frequencies where also most of the signal is lo- 
cated. This implies that it will be impossible to perfectly sep- 
arate the true signal from the large-scale structure contamina- 
tions on the same scales. On the other hand, the lensing noise 
contributed by the large-scale structures on scales larger than 
the typical scale where the signal power spectrum becomes 
negligible will be effectively filtered away. 

4. Aperture mass 

For detecting halos throu gh their weak-len sing signal, the aper- 
ture mass statistic (M^p. ISchneideil [19961) is frequently used. 
Due to its favourable statistical properties, it has turned into 
one of the prime tools for analysing cosmic-shear statistics. Its 
application to the detection of the weak-lensing patterns from 
dark-m atter halos has been shown by several authors (see e.g 
lErben et al.. 2000; Schiriner et al, ,2004.) . 

The aperture mass is a weighted integral of the local surface 
mass density, 

Map(6/) = J K{0')U{e-ff) , (33) 

where U{0) is the weight function. If the aperture is chosen to 
be circular, U is axially symmetric, U{0) - U{6). Since the 
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convergence is not directly measurable, Eq.[321is conveniently 
rewritten in terms of the tangential shear yt as 



Schneider 1996 - 

Optimal filter, excluding LSS - 

Optimal filter, including LSS r^^2 arcmin ■ 

Optimal filter, including LSS fg = 4 arcmin . 



provided the filter U is compensated, i.e. satisfies 

/ 

Jo 



d0' 0' U(0') = , 



The function Q{0) is related to U{0) by 



Q(0) 



" 0' Jo 



de' 0' u(0') - u(0) . 



The noise of the aperture mass is given by its variance 



ng Jo 



A0'0'Q^{0') 



(34) 



(35) 



(36) 



(37) 



llSchneidei{ll996l) . 

The function U can be chosen such to maximise the sienal- 
to-noise ratio M^^Ictm.^^- Schneider et al. (1998) show that this 
is the case if Q mimics the tangenti al shear profile of the lens. 
For example, ISchrrmer et a n ll2004 use a function which ap- 
proximates the expected shear profile of a NFW halo. 

Thus, an assumption on the expected profile of the tangen- 
tial shear needs to be made in the construction of both the op- 
timal filter and the aperture mass. In this respect, the two func- 
tions U and T should be chosen to look very similar. Two differ- 
ences are important, however. First, the weight function U used 
in the definition of the aperture mass is compensated, while t 
does not need to be. The reason is that the aperture mass was 
specifically defined such as to be proportional to the projected 
mass enclosed by the aperture. Our goal here is more modest 
as we only ask for significant halo detections, rather than for a 
determination of halo parameters apatr from the overall shear 
amplitude. This allows us more freedom in the choice of t com- 
pared to U. Second, the shape of the optimal filter differs from 
the profile function t because it is modified by the shape of the 
total noise power spectrum. A comparison between the optimal 
filter and the function Q{0) used by Schneider (19_9^ is shown 
in Fig. 13 The solid and the dot-dashed lines refer both to the 
optimal filter but for two different scale radii used for mod- 
elling the signal. We will show in the following sections how 
the two filters perform when apphed to simulated data. 

5. Simulated source catalogues 

Numerical simulations are a unique tool for assessing the per- 
formance of the optimal filter, individually and in comparison 
to the aperture mass. In the following sections we describe the 
lens models and the numerical methods used for producing 
simulated catalogues of weakly distorted sources, subsequently 
analysed with our filter and the aperture-mass statistic. 

5.1. Numerical models 

The numerical model used here is a simulation of a super- 
cluster hke region containing six halos with virial mass larger 




Fig. 2. Comparison between the filter profiles. Different lines 
refer to the conventional aperture mass (solid line), the optimal 
filter described in this paper but without the large-scale struc- 
ture contribution to the noise power spectrum (dashed line), 
the full optimal filter for r., = 2' and = 4' (heavy solid heavy 
dot-dashed lines, respectively). The large-scale structure power 
spectrum modifies the filter shape according to the filter scale, 
by lowering its sensitivity on large scales where weak lensing 
by large-scale structure dominates. 



than 5. 6x10'^/;' Mq. They are embedded in an over-dense re- 
gion containing several less massive halos. The region is a sec- 
tion from an existing dark-matter only large-scale cosmologi- 
cal simulation, later re-simulated at higher resolution usin g the 
"Zoomed Initial Conditions" (ZIC) technique ( Tormen et all 
[l997). The parent simulation was an A^-body run of a box 
of 479 /i ' Mpc side length containing 512^ particles of mass 
6.8x10'" /i"' Mq. The background cosmological model is spa- 
tially flat with fim = 0.3 a nd Qa = 0.7 at the final epoch, iden- 
tified with redshift zero jYoshida et all l200l[ I.Tenkins et all 
.2001). 

From the output of this large cosmological simulation at 
redshift zero, we selected the super-cluster region for the fol- 
lowing re-simulation. New initial conditions were constructed 
as follows. First, the particles contained in the selected region 
were traced back to their initial (Lagrangian) positions. These 
define a region in Lagrangian space. Then, the initial density 
field in that Lagrangian region was re-sampled by placing a 
higher number of particles than were originally present, adding 
additional small-scale power appropriately. Gas was introduced 
into the high-resolution region by splitting each parent parti- 
cle into a gas and a dark-matter particle. The final mass res- 
olution of these simulation was otdm = 1.13 x 10^/i"'Mq and 
'Wgas - 1 .7 x 1 0'^ ' Mq for dark matter and gas within the high- 
resolution region, respectively. Thereby the simulation hosts a 
usable region of 25 Mpc. 

Finally, the region was evolved down to redshift z - us- 
ing the code GADGET-2, a new version o f the parallel tree-SPH 
simulation code GADGET d Springel et aL. 2001) . The simula- 
tion was carried out including only non-radiative hydrodynam- 
ics. The comoving gravitational softening length was fixed to 
e = 5.0 /i ' kpc (Plummer-equivalent). 
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For the present work, we use the output of the re-simulation 
at redshift z - 0.297. We separately analyse six sub-boxes, 
each of them centred on one of the six most massive ha- 
los. The comoving side-length of each sub-box was set to 
~ 14.45 /j"' Mpc, corresponding to 1° on the observer's sky. 

5.2. Weak-lensing simulations 

The particles in each sub-box sample a three-dimensional 
density field. Their positions are interpolated on a grid 
of 512^ cells using the Tria ngular Shaped Cloud method 
llHocknev & Eastwoo 1 Il988h . Then, the three-dimensional 
density field is projected along its coordinate axes, thus pro- 
ducing three surface density maps, which we use as indepen- 
dent lens planes in the following lensing simulations. 

The lensing simulations are performed by tracing a bundle 
of 2048 X 2048 light rays through a regular grid, covering the 
central quarter of the lens plane, corresponding to a field of 
view of 1°. The deflection angles a^igii, j) of each ray (/, j) are 
computed acco rding to the method described in sev eral earlier 
papers (see e.g. lMeneghetti et allEoO O. 2001, 2004). 

Weak gravitational lensing by large-scale structures be- 
tween the observer and the sources is included as follows. From 
the power spectrum i26i of the effective convergence, we ob- 
tain the power spectrum of the effective lensing potential. 



P^(k) = -P,(k) . 



(38) 



Since we only aim at measuring the efficiency of our filter for 
detecting a shear signal, the source redshift distribution is not 
important. Thus, all sources are assumed to be at the same red- 
shift Zs, considering the two cases Zs = 1 or = 2. The upper 
integration limit in Eq. j23> is then the comoving distance to 
the source plane, - w(zs), and the distance distribution of 
the sources is a delta function. 



G(w) = 6(w - Ws) 



(39) 



We generate the effective lensing potential in Fourier space 
as a Gaussian random field. Using Eq. (|2j, we find the large- 
scale structure contributions to the deflection angles, QrLSs(', j)- 
The total deflection angle for each ray (/, j) is then 



fftot = ffsig + ffLSS 



(40) 



Reduced-shear maps for each field are straightforwardly de- 
rived from the deflection angle as explained in Sect.|2 

Background galaxies are randomly placed and oriented, 
and their ellipticities drawn from the distribution 



Ps(ksl) = 



exp[(l-hp)/^2j 
;rcr|[exp(l/(r|)-ll 



(41) 



with cr^^ = 0.3. We assume a background galaxy number 
density of rig - 25 arcmin"^. Therefore, each field contains 
~ 90, 000 galaxies. "Observed" ellipticities are obtained from 
Eq. il 1> . and a catalogue containing the shapes and the posi- 
tions of all sources in each field is compiled for the following 
analysis. 



6. Results 

6.1. Applying the filter 

We now describe how we analyse the simulated catalogues. 
The integral in Eq. (I15> is approximated by a sum over galaxy 
images. Since the ellipticity e is an estimator for ly, we can 
write 



Aest(0) = - Ve^m^vm - 0\) 



(42) 



where euiOi) denotes the tangential component of the /th 
galaxy's ellipticity relative to the line connecting 0, and 0. 
Thus, the filtering is performed in the real domain, which 
requires to Fourier back-transform the filter function from 
Eq. (HIJ. 

We measure Aest on a regular grid of 128 x 128 points cov- 
ering the field. The noise estimate in Aest is obtained from 



cr 



(43) 



We produce maps of the signal-to-noise ratio (S /N) for all our 
fields. We then search for peaks with heights above a minimal 
S/N. 

For comparison, similar signal-to-noise maps are produced 
using the aperture mass statistic. In this case, we also use 
Eqs. J42> and J43> for estimating the signal and the noise, but 
we substitute the filter function Y with the function Q discussed 
in Sect.m 

6.2. Signal-to-noise maps 

We show in this section how the filter performs on one projec- 
tion of the field centred on the third most massive halo in the 
simulation box. The results for the fields centred on the other 
five halos are qualitatively and quantitatively very similar. The 
simulation we use represents a fairly complex case to study. 
Indeed, the field is crowded with ~ 15 subhalos of different 
masses and some filamentary structure connecting them. 

Both the optimal filter and the aperture mass are designed 
for the optimal detection of halos of a given scale. For the opti- 
mal filter, this can be set to the scale radius of the NEW profile. 
However, this definition does not apply to the aperture mass. 
In fact, the size of the aperture is defined as the radius where 
the filter function Q has its first root. In order to compare the 
efficiency of the two filters for the detection of the same class 
of object, we have empirically adjusted the size of the aperture 
such that it resembles the scale of the signal profile used in the 
optimal filter. 

We show in Fig. |3 the signal-to-noise contours for detec- 
tions reaching above 5<t. The left and the right panels refer to 
the aperture mass and the optimal filter, respectively. From top 
to bottom, the filter scale is 1', 2', 4' for the optimal filter, corre- 
sponding to scales of 2.75', 5.5' and 11' for the aperture mass. 
In order to allow the reader to see which detections are real 
and which are spurious, we overlay the contours on surface- 
density maps from the numerical simulation. The halos which 
are identified by the halo finder are marked by circles with radii 
proportional to their virial radii. 
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Fig. 3. Signal-to-noise contours above 5cr for different filter 
scales, assuming sources at redshift Zs = 2. Results obtained 
by applying the optimal filter and the aperture mass are shown 
in the left and right panels, respectively. From top to bottom, 
the scales are 1', 2' and 4' for the optimal filter, and 2.75', 5.5' 
and 11' for the aperture mass. The contours are overlaid in all 
panels on the surface-density map in the field centered on the 
third most massive halo in our simulation. 
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Fig. 4. Signal-to-noise contours for a fixed filter scale of 2' for 
the optimal filter (left panels) and of 5.5' for the aperture mass 
(right panels) for sources at redshift Zs - 2. From the top to 
the bottom panel, contours start at 2cr, 4cr and 6cr, respectively. 
In all panels, the contours are overlaid on the surface-density 
map in the field centered on the third most massive halo in our 
simulation. 



The optimal filter performs much better than the aperture 
mass in detecting the dark matter concentrations present in 
our simulated field. On all scales, the signal-to-noise ratios are 
larger for the optimal filter than for the aperture mass. 

Spurious detections occur with both filters. A further in- 
spection shows that these false matter concentrations are in fact 
connected to over-densities produced by the projection of the 
large-scale structures along the line-of-sight. As discussed in 
Sect. 13.31 these misidentifications cannot strictly be avoided, 
for the large-scale structure power spectrum is still large on 
the spatial scales where most of the signal resides. A possible 
method for discriminating between true halos and large-scale 
structure superpositions could be to group the sources in red- 
shift bins. In fact, the large-scale structure pattern is expected 
to change for sources at different redshifts. This method is cur- 



rently being tested using numerical simulations (Meneghetti et 
al., 2005, in preparation). 

It is remarkable that large-scale structure contaminations 
are much less important for the optimal filter compared to the 
aperture mass, especially when aiming at detecting clusters of 
large mass and low redshift, having large scale radii. 

In Fig. 13 we show the detections above a minimal signal- 
to-noise ratio for a fixed filter scale of 2 arc minutes for our 
filter (right panels) and of 5.5 arc minutes for the aperture mass 
(left column). From the top to the bottom panels, the contours 
start at 2, 4 and 6cr, respectively. Again, the efficiency of the op- 
timal filter is remarkably higher than that of the aperture mass. 
No detections above 3cr are found with Map, while several true 
peaks are found using the optimal filter at more than 6cr. 

Apparently the chance of a halo to be detected does not de- 
pend only on its mass or on the amplitude of the shear which 
it is able to produce. Indeed, the noise produced by the large- 



8 



M. Maturi et al.: An optimal filter for the detection of galaxy clusters through weak lensing 



scale structure can locally disturb the signal, such that even a 
relatively large halo can fall below detectability due to the local 
effective shear produced by the large-scale structure. As a sim- 
ple test, we perform several weak-lensing simulations keeping 
the background galaxy population fixed and changing the large- 
scale structure in front of them, by using different realisations 
of the Gaussian random deflection-angle field. The results are 
shown in the sequence in Fig.|5]for sources redshift Zs - 2, and 
in Fig.|6lfor Zs - ^- We use a filter scale of 2 arc minutes and 
plot the signal-to-noise contours above 4cr in all panels. In the 
right panels of images, the optimal filter is applied in order to 
remove the large-scale structure contamination. Clearly, some 
of the clumps which are detected for a given realization of the 
large-scale structure, are entirely missed when a different seed 
is used. For better illustrating the influence of the large-scale 
structures on the detection efficiency, we show in the left panels 
the same detections without including the large-scale structure 
noise in the construction of the optimal filter. Spurious halo de- 
tections are strongly reduced by using low-redshift sources be- 
cause of the weaker large-scale structure contribution along the 
line-of-sight. This causes the signal-to-noise ratio to increase 
despite the clusters' being less efficient lenses for the lower- 
redshift sources. Compared to Fig.|5]the clusters' lensing effi- 
cency is smaller because of the lower sources distance, but, for 
the same reason, there is a strong reduction in the large-scale 
structure contamination which leads to a gain in the signal-to- 
noise ratio. 

7. Summary and discussion 

We have developed a linear filter for detecting the weak-lensing 
signal of dark-matter halos. Given the expected shape t of the 
halo signal, the filter is uniquely determined by the two require- 
ments that it be unbiased and maximise the signal-to-noise ra- 
tio. The filter shape is modified by the noise power spectrum, 
into which we include the shot noise from the random posi- 
tions of the galaxies used to measure the gravitational shear, 
the noise from their random intrinsic ellipticities, and the weak- 
lensing signal contributed by the large-scale structure. Since 
the total noise power spectrum is monotonically falling over 
the scale of the shape function t, the inclusion of the noise nar- 
rows the filter considerably. 

Our primary goals are to obtain a filter which separates ha- 
los from spurious peaks in the lensing effects of the projected 
large-scale structure, and whose results are stable against small 
changes in the filter scale. Our motivation is that applications 
of another statistic used for detecting dark-matter halos, the 
aperture mass Map, tend to detect large numbers of dark peaks, 
i.e. peaks in the lensing signal which are not found to corre- 
spond to detectable signal in the optical or X-ray wave bands. 
Moreover, these peak detections tend to be unstable against 
changes in the characteristic scale of Map, i.e. the significance 
of their detection may change wildly upon even small changes 
of scale. 

It is clear that a clean separation of halos from large-scale 
structure is impossible because large-scale structure can be 
considered as being composed of dark-matter halos. We thus 
choose the operational definition that the halos we wish to de- 
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Fig. 5. Signal-to-noise contours above 4cr for different realiza- 
tions of the large-scale structure. Results obtained by applying 
the optimal filter with and without filtering out the large-scale 
structure are shown in the left and in the right panels, respec- 
tively. The scale of the filter is 2' and the source redshift is 
Zs - 2. 

tect are distinguished from large-scale structure as being ob- 
jects on scales smaller than the non-linear scale. Consequently, 
we model the large-scale structure "noise" in the weak-lensing 
signal as being due to linearly evolved structures alone. 

We perform numerical simulations for testing our new fil- 
ter. The results, summarised in Figs. |3] E] and |5] demonstrate 
that the filter meets all of our goals. Changes in the filter scale 
leave the signal much more stable than the aperture mass does, 
as Fig. |3] shows. Fig. |4] illustrates that the signal-to-noise ratio 
produced with our filter is typically substantially higher than 
that of the aperture mass, and Fig.|5ldemonstrates that filtering 
halos against large-scale structure lensing is indeed very effi- 
cient. Compared to the aperture mass, the number of spurious 
peaks is also substantially reduced. 

We thus propose this filter as an alternative statistic for de- 
tecting dark-matter halos. Of course, the aperture mass remains 
a very valuable statistic which has its significant strengths in 
analyses of lensing by large-scale structures. It appears, how- 
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Fig. 6. As Fig. Is] but with the sources placed at redshift Zs = 1 ■ 

ever, from our numerical simulations that a fair fraction of the 
peaks found with the aperture mass may in reality not be caused 
by halos, but by peaks in the projected large-scale structure. 
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